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Abstract 

The Selberg correlation integrals are averages of the products Y\™=i YYl=i( x s ~ z i)^ a with 
respect to the Selberg density. Our interest is in the case m = 1, m = fj,, when this 
corresponds to the /x-th moment of the corresponding characteristic polynomial. We give 
the explicit form of a (n + 1) x (n + 1) matrix linear difference system in the variable 
which determines the average, and we give the Gauss decomposition of the corresponding 
(n + 1) X {n + 1) matrix. For a positive integer the difference system can be used to 
efficiently compute the power series defined by this average. 



1 Introduction 

The probability density function (PDF) on z\ G [0, 1] (/ = 1, . . . , N) 

1 n 

— ^n^a-^r- 1 n i^-^i 2t > 



S , „(ai,a 2 ,r) 

1=1 iSj<k<n 



where 



„ n 

S n ( ai ,a 2 ,r) := / JJ - z^ 1 J] \z 3 - z k \ 2r dz x ■ ■ ■ dz n 

JlOA}™ i=l l<j<k<n 

TT r(q 1 +jr)r(a 2 +jr)r(l + (j + l)r) 

f\ r( ai + a 2 + (n + j-l)r)r(l + r) ' [) 

is the Selberg integral, plays a fundamental role in the theories of random matrices and 
Calogero-Sutherland quantum many body systems (see e.g. [8]). In random matrix theory, 
(1) with (a 1 ,a 2 , r ) = (a + l,b + 1,(3/2) defines the Jacobi /3-ensemble. For (3 = 1,2 and 
4, and certain a, b, this can be realized as the eigenvalue PDF occurring in the analysis of 
correlation coefficients associated with Gaussian data sets, and also as the singular values of 
sub-matrices formed from various unitary matrices chosen according to the Haar measure. For 
general ai,a 2 ,T > there are constructions of (1) relating to similarity reductions of unitary 
matrices to Hessenberg form [17], to block diagonal form [5], and to the generalized eigenvalue 
problem for certain tridiagonal matrices [10, 24]. 
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Upon the change of variables Zi = sin 2 0j, < fa < ir/2 (i = l,...,n), the PDF (1) 
becomes proportional to 

n 

J](sin 2 fap (cos 2 fap J] | sin 2 fa - sin 2 fa\ 2r 

i=l l<j<k<n 

with a[ — ai — 1/2; oc' 2 = a 2 — 1/2. As such it is the absolute value squared of the ground 
state wave function for the .BC-type Calogero-Sutherland Schrodinger operator [8, Eq. (11.55)] 

_ _|_ f a i T ( a i T ~~ x ) + a 2 T ( a 2 T - 1' 



+2r(r- 1) 



sin' 0i cos 2 < 



i<j^k<n Vsin2 (^ - <t>k) sin 2 ^- + fa 



It is well known that (1) exhibits many remarkable integrability properties. One is the 
gamma function form of the normalization (2). Another is that the family of averages associ- 
ated with (1) 



/ / flz^il-Zif-^x-z^ TT \z j -z k \ 2T dz 1 ---dz n 

J[0,z]> J [0,1]"-* i=1 Ki<k<n 



(3) 

l<i<fe<" 

(we have set a.\ = a, a 2 = (3) can be characterized in terms of a certain differential-difference 
equation [7], equivalent to a (n + 1) x {n + 1) matrix Fuchsian differential equation [11, 20]. 
With p = and fi = 2r, (3) is simply related to the one-point density implied by (1), and the 
differential-difference equation was used in [7] to compute the polynomial in x specified by (3) 
in the case r G Z + (the polynomial is of degree 2rn so for practical purposes rn cannot be 
too large). 

In the case r = 1 (3) can be calculated in terms of the solution of the Painleve VI non-linear 
differential equation in a form [12]. It is also revealed in [12] that the a-function associated 
with (3) satisfies a third order non-linear difference equation for integer shifts in the variable 
/i, while (3) itself can be computed by a recurrence scheme based on the discrete Painleve V 
equation. This can be understood from the viewpoint of a more general theory relating to 
isomonodromic deformations of linear difference systems [3]. 

It is the objective of this work to provide a (n + 1) x (n + 1) matrix linear difference system 
for integer shifts of the variables a, (3 or // (the latter restricted to cases that (x — Zi)\x — = 
±\x — for some sign ±, or alternatively to twice integer shifts) in the integrals (3) with 

r > 0. This will be used to provide an alternative method to compute the polynomial in x 
specified by (3) in the case p — 0, /i — 1 even. 

A precise formulation of the family of Selberg correlation integrals to be studied is given in 
Section 2, along with a statement of our result for the explicit form of the difference system. In 
Section 3 we introduce a certain family of interpolation polynomials, and we state three term 
relations satisfied by the corresponding generalized Selberg integrals. The three term relations 
are shown to imply the difference system. We give their proof in Section 4. In Section 5 it is 
shown how to use the difference system to compute (3) in the polynomial case. Furthermore, 
we specify applied studies in random matrix theory to which our computations have relevance, 
and we make note too of the wider problem of characterizing correlation functions in statistical 
mechanical problems in terms of differential or difference equations. 
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2 Definitions and main result 

We begin with some definitions. We introduce the notation $(z), z := (z±, . . . , z n ), to denote 
the generalization of the integrand (3), 

n 

*(z) := J] \ Xl - ^| ai_1 k2 " z^lxs - z^- 1 J] \ Zj - z k \ 2T . (4) 

i=l l<j<k<n 

The parameters a±, a 2 , 013 are assumed restricted to domains for which it is possible to specify 
a region A C M. n with the property $(z) vanishes on the boundary dA of A. For example, if 
< x < 1 and Re(ai), Re(a2), Re(a3) > 0, we can specify 

A = A p = [0 l i]'x[0,ir (p = 0,l,...,n). (5) 
For rational functions 0(z) bounded on A we define 

(0) := J <j>(z)$(z)dz 1 ---dz n , (6) 

and we use this notation in turn to specify T a . according to 

n 



1=1 



Note that in the cases that — Xj)\zi — Xj\ ai_1 = ±\zi — Xj\ aj for some sign ±, T aj corresponds 
to incrementing otj by 1, and independent of this requirement, T% corresponds to incrementing 
aj by 2. Our goal is to identify polynomials {<Pi(z)} i= o,i,..., n +i such that {T ai }i=o,i,...,n+i is 
linearly related to {(</?*) }j=o,i,...,n+i- 

For this purpose we take inspiration from the work of Aomoto [1, 2]. With <S>*(z) denoting 
®(z) specialized to x\ — 0, x 2 — 1, ql% — 1 so that 

n 

^w^n^a-^r- 1 n i^-^i 2t (8) 

i=l l<i<fc<n 



and with 

(0)*:= / (9) 

7 [0,1]" 

it was proved in [2] that 



i+l 



ai + «2 + (2n — % — 2)t 



Since S n {a.\ + 1, a 2 , r) = § z\ - • ■ z n $*{z) dz\ • ■ ■ dz n , by iterating this we immediately have the 
difference equation for the Selberg integral of (2), 

+ l,a 2 ,r) = S(a 1 ,a 2 ,r)f\ Ql ~j~ ^ — ^ , (10) 

±J - ai + «2 + (2n - i - l)r 

«=i v 7 
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which in turn can be used to deduce the Gamma function evaluation given in (2). Thus we 
learn that in the case of the Selberg integral T ai (l}* is linearly related to (1)*, and moreover 
we note that to derive the linear relation use was made of the auxilary polynomials {111=1 z i} 
- referred to as interpolation polynomials for the role they play in the calculation. 

We will follow the same general strategy in relation to the integrals (6). Thus a family 
of interpolation polynomials, {<fi,j(z)}, will be introduced so as to determine the polynomials 
{</?j(z)}j=o,i,...,n+i forming a polynomial basis of the difference system associated with the shift 
(X\ i — y (X\ ~\~ 1 in the integrals (6). We know from [1, 2] that the main tool in determining these 
polynomials is the vanishing of a certain class of averages (6). 

Lemma 2.1 For k = 1,2, ... ,n let 



— 1 chq — 1 v— v 2r \ . . . 

— Zk X 3 — Z k £j£ n Z k — Z[/ 

We have (V fe 0) = 0. 



dz k <b(z) dz k 

Oil — 1 

X\ — z k 



Proof. By definition, 

(V fe 0) = J <5>(z)V k (i)(z)dz 1 ---dz n = J ■^-(^>(z)^(zfjdz 1 ---dz n = 
if (f)(z)§(z) vanishes on the boundary dA of A, which we have assumed. □ 



However, for purposes of presentation, rather than to start with the interpolation polynomi- 
als, it is convenient to immediately present our findings for the explicit form of the polynomials 
{</?i(2)}i=o,i,...,n+i and the corresponding difference system. 

Theorem 2.2 Write 



ifii(z) := (x 2 - z ± ) ■ ■ ■ (x 2 - z n -i) x (x 3 - z n _ i+1 ) ■ ■ ■ (x 3 - z n ) i = 0, 1, . . . , n. (12) 



— v — 

' 1 



We have 



T ai ({(p ), {(fl), ((p n )) = ({(p ), (Vl), • • • , (<Pn))A 



where A = LDU with 

/loo \ 
ho hi 

\ hiO Inl ' ' ' Inn ) 



/d 



D 



d 1 



\ 



d n J 



, U 



(13) 



/ -u 00 ""oi • 

Mil ■ 

V 



U 0n \ 
Uln 

Unn ) 
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All entries in L, D, U not explicitly shown are zero, while for the non-zero entries we have 



I 



-1 



(a 2 + jt; t 



i-3 



X 2 - Xi 



n-j 

n-ij (ai + a 2 + 2jr; t)_\x 3 - x x 

' 1 3 



(ax\ r)j(o;i + a 2 + 2jr; r) n _.(x 2 - Xi) J (x 3 - x x ) n 1 



(ai + a 2 + (j - l)r; r) (ai + a 2 + «3 + (« + j - 1)t; r) 



I/, 



j\ («3 + (n - j)r; r) 



where (x; r) = 1 and (x; r)j := x(x + r)(x + 2r) • • • (x + (i — l)r) /or i = 1,2, . . .. 



(14) 



The proof of this result will be given in Section 3. 

In Theorem 2.2 the matrix A is given in terms of its Gauss LU decomposition. The 
symmetry of (4) under the interchange (x 2 ,a 2 ) i— > (x 3 ,a 3 ) allows for A also to be written in 
terms of its UL Gauss decomposition. To see this, let L, D and U be the matrices L, D and 
U after this interchange. Since 



«¥>n), (Vn-l), • • • , (v?o)) = ((<A)), (^l), • • • , 



where 



J 



/ 



V 1 



1 \ 



/ 



(4) can be rewritten 

T ai ((<p ), M, (<p n )) = ((<p ), (if,), (<p n ))U'D'L', 

where 

U' = JLJ, D' = JDJ, U = JUJ. 

We see that U' , D' and V are upper triangular, diagonal and lower triangular matrices, 
respectively with non-zero entries 



< 3 = i \y :( 



(«3 + (n - j)r; r 



(ai + a 3 + 2(n - j)r; r) - a>i 



x 3 - Xi 



c£ = 



I- 



(ai, r) n _j(oi!i + a 3 + 2(n - j) r ; r^fo - zi) J (x 3 - ^i) n J 
(ai + a 3 + (n - j - 1)t; ^) n _ i («i + a 2 + «3 + (2n - j - l)r; r)^. 

(a 2 + jr;r)._. 



- V («i + «3 + 2(n - i)r; r) 
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3 Interpolation polynomials 



In this section we present some lemmas, and collaries of the lemmas, which together imply 
Theorem 2.2. We will defer the proof of one of these — certain key three-term relations - 
until the next section. We begin with a lemma which enables the difference system of Theorem 
2.2 to be written in a more convenient form. 

Lemma 3.1 Let U be as in Theorem 2.2. We have that U^ 1 = (u*j)o<i,j< n the upper 
triangular matrix with non-zero entries 



u 



Ki 



(a 3 + (n- j)r; r) 



i) [a l + a 2 + {] +i-l)r;r) j _.' 
Proof. It suffices to check that for i < j, J2k=i u ik u tj = 0- Now 



Yl UikU *kj 



k=i 



t_i fk\ ("3 + (n - k)r, r) t _. / j\ (a, + (n - j)r\ t) 



k=i 



j-k 



(«! + a 2 + 2ir; r) \kj (on + a 2 + (j + k - l)r; r) 



j-k 



(«3 + (n-j)r;r) 



(«i + a 2 + 2ir; r) 



2j-2i-l 
j-i 



(ai + a 2 + 2ir; r) 



2j-2j-l 



k=i 



(a 3 + (n-j)r;r)._. fJ ,^_^ k(J _ l 



k=0 



. ,(ai + a 2 + 2ir + (k — i)r; t) . . . 
k — %< 

(a.\ + a 2 + 2ir + kr; r) , 

V / 1 % 1 



fc 



and the last summation vanishes as an example of the summation formula for 2-^1(0, b; c; 1). □ 
A more convenient form of the difference system can now be established. 



Lemma 3.2 



( Uoo Uqi 
Mil 

V 



U n\ 



((<^o),(^l),---,(^n)) 



/ [00 

^10 hi 



Inn ) 



where 



u 



(«i + a 2 + a 3 + (n + j - l)r; r) (a 3 + (n - j)r; r) («i + a 2 + (j - l)r; r) ., 



(-1) 



t-3 



n-j 



n — i 



[on; r)j(a 2 + jr; r)i_j («i + a 2 + (i + j) r ; ^^(a* - ^i) n *(^2 - xi) 1 



This is obtained by acting on the left of both sides of the difference system of Theorem 2.2 
by C/ _1 , making use of the explicit form of the latter known from Lemma 3.2 on the LHS, and 
clearing denominators. 
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Before considering the proof of this rewrite of the difference system, we make note that in 
a special case it implies the recurrence (10) for the Selberg integral. Thus we note from (12) 
that with x 2 = %3 we have <Pi(z) = YYi=i( x 2 ~ z i) independent of i. We note too that it follows 



from the definitions of and 1^ in Lemma 3.2 that 



3 



^ Uij = u o = («i + oc 2 + a 3 + (n - l)r; r)„ 

i=0 
n 

^ kj = Ln = (Cti; r) n (x 3 - X 1 ) U , 



1=3 



\X2=X\i 



valid for j — 0, 1, . . . , n, where use has been made of the summation formula for 2-^1(0, b; c; 1). 
It follows that if we set x 1 = 0, x 2 = x 3 = 1 and replace a 2 + a 3 with a 2 in the difference 
system of Lemma 3.2, then it degenerates to a single equation, which is precisely (10). 

To derive the difference system of Lemma 3.2, and thus that of Theorem 2.2, we introduce 
the interpolation polynomials <Pi,j(z) according to 

<Pij(z):=(z 1 -x 1 )(z2-x 1 )~-(z j -xi)(pi(z) for i,j = 0, 1, . . . , n. (15) 

V v ' 

3 

Note that setting j = gives (fi,o(x) = V 9 ^- 2 )) while setting j = n gives 

<Pi,n{ z ) = ( z i - x i)( z 2 -xi)--- (z n - xi)(pi(z), i = 0, 1, . . . ,n 

and so 

T ai (<Pi) = (<P*,n) , (<P%) = {<Pi,o)- (16) 

Most importantly, the integrals (6) with = (pi+ij, <Pi,j+i, <fi+i,j+i, or with = (pij+i, <fij, 
(fii+ij satisfy certain three-term relations. 

Lemma 3.3 (Three-term relations) For i,j — 0, 1, . . . , n — 1 we have 

(ax + (n-j- 1)t)(x 2 - x 1 )(^ +lj ) 
= (a 3 + {n - i - 1)t) (<p itj+1 ) + («i + a 2 + (n + i - j - 1)t) {(p i+1>j+1 ) . (17) 

(ai + a 2 + a 3 + (2n - j - 2)r) (v^ij+i) 
= (ai + a 2 + (n + i - j - l)r)(x 3 - ^i)(^j) - (a 2 + ir)(x 2 - xi){(p i+1J ). (18) 

Proof. See Appendix A. □ 

These three-term relations in fact imply the difference system of recurrences in Lemma 
3.2. To see this we first use an induction on j to deduce from the three-term relations two 
particular difference systems. 
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Corollary 3.4 For 0<j<k<nwe have 

(an + (k - j)r;r) .(x 2 - x^ 3 ((p k ^ k ) (19) 

= { ) (" 3 + ( n ~ k "> T > T )j-i( a i + a 2 + (2fc - j - l)r;r).(^ + fe^ >t _fc +i ), 

i=o 

(«i + a 2 + «s + (n + j - 1)t; r) n _ jfc (^ >n _ J ) (20) 

= HC" 1 )* - * [ _ • («i + «2 + (? + 2j - fc)r; r) n _.(a 2 + jr; r),_ fe 
i=fc Vn 

X(X 3 - X l ) n ~ i (x 2 - XiY^i^i-k+j^-j). 

In particular, by setting k — j in the above, for j — 0,1, ... ,n we have 

(ai;r)j(x2-a:i) J <¥>j,n-;j) (21) 
= l J J ("3 + (n- jO^r)^.^! + a 2 + (j - l)r;r).(v?i,„), 

(«i + a 2 + a 3 + (n + j - l)r; r) ^((fij^-j) (22) 
= E(-!) i_J (" ~_[) ( ai + "2 + (» + j)r; r) n _.(a 2 + jr; r)< 



x(x 3 - Xi) n *(x 2 - an)* j (^,o>- 

Proof. See Appendix B. □ 

The difference system of Lemma 3.2 can now be derived. 
Proof of Lemma 3.2. Multiplying (21) and (22) by appropriate factors so as to make 
their LHS's equal, then equating RHS's gives 

j n 

^Uiji^n) = ^lij&ifl), j = 0, 1, . . .,71, 

i=Q i=j 

where Uij and Uj are as in Lemma 3.2. Making use of (16) this reads 



i=3 



which is precisely the sought difference system. □ 

4 Implementing the recurrences 

Consider $ specialized to the function <J>*(z) of (8) but with a\ = a, a 2 = (5. Let this be a 
consequence of setting 

cki = 1, a 2 = a, a 3 = (3, x± — x, x 2 = 0, x 3 — 1 (23) 
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in (4). With this choice of <£>* define = ((f))* /(l)*, where (</))* is specified by (9). Here 
our aim is to use the difference system of Theorem 2.2 to explictly compute the polynomial in 
x specified by 

Let (T ai ((p })# refer to (7) with the substitutions (23) made afterwards, and normalized by 
dividing by (1)*. Then according to (7), (12) 

^w^^^((n^-^) # L„J- <»> 

Thus our task is to compute (T^^q))*, as we have 

(Wf = { _ ir ^s_^M (( _ , r )*|_J. (26) 

For /i = 1 the closed form evaluation of (24) is known from the work of Aomoto [2] as being 
proportional to the Jacobi polynomial P^ /1 '"' 2 \l — 2x) with 7x = a/r — 1, 72 = /3/t — 1. This 
in turn can be written in terms of a Gauss hypergeometric function, giving 

/ n \* 

(J[(x-Zj)) = c 2 F 1 (-n,(a + p)/T + n- l,a/r;x), (27) 
j'=i 

where 

g = (-m^r) n 

(a + /3 + (n-l)r;r) n V ; 

(the factor of c is required to make the coefficient of x n on the RHS unity). We can use 
knowledge of this to calculate (y?fc) # (k — 0, . . . , n). Once these have been determined we can 
use (13) with A =: A ai to recursively compute (T^^o))* according to 

(t£«^ >, W, • • • , (V^n)) # = ((^0) # , • • • , (</>n> # Ml^2 • • • A„ (29) 

Upon making use of (26) this determines (24). Explicitly, with (v)k denoting the fc-th compo- 
nent of the row vector v, we have 

(- i )" < "" 1) to^t « jWD LJ - (<«*. ■ ■ ■ ■ A "){ 

(30) 

Lemma 4.1 Le£ fre specified as below (23). We have 

, N# , lVt (a;r) w (-/3-(n-l)r;r) fc 

(V?fc) = (a + /> + (n-l)r;r) B ' (31) 
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Proof. According to (7), (12) and (27) 

(T n <V?o» # = gw ^ +1 /\ r) c 2 F 1 ( - n, 1 - fj/r - n, a/r; --?-) , (32) 

where use has been made of a Kummer relation for 2 Pi- Note from (10) that 

S n (a + l,/3,r) (a; r) n 

S n (a,(3,r) (a + /3 + (n- l)r) n ' 

On the other hand, it follows from (13) that 

n 

(T ai {p )f = d J2&k)*k (33) 

fc=0 

where, after substituting (23) in (14), 

(1 + a;T) n (a;r) k 



d l 



kO 



(1 + a + (5 + (n - 1)t; r) n (l + a; r) k \n - k) VI - x) 

Equating (32) and (33) we see that the factor of (1 — x) n cancels, and we can equate coefficients 
of (-x/(l - x)) k to deduce (31). □ 

It is also possible to derive (31) independent of knowledge of (27), using instead an in- 
tegration formula in Jack polynomials theory, due to Warnaar [21] (see also [6]). With 
A = (Ai,...,A„) a partition of non-negative integers, and P^ a \t) = pjf*\ti, . . . , t n ) denot- 
ing the symmetric Jack polynomial, the integration formula reads 

„ n 

/ pi i/r) wn^r i (i+^)" a: " 2/ " 2(n " i)T n \t k -t j \ 2r dt 1 ---dt n 

J[0,oo) n i= l l<j<k<n 



P? M W\ '^'"-''^ - S n{x ,y,r), (35) 

\ Vl =-=V n =-l \ — y _|_ 1 1 \ I > 



= ...=„„ = -! \_y _|_ ]_] 



(1/t) 
A 



where 



[< r(«-a-i)/a) • (36) 

Our interest is in a transformed version of (35). 

Corollary 4.2 Let (<f>)# be as specified below (23) and letP^ /T \~^) = pj 1/r) (^, . . . , 



[0 + (n - 1)t]< 1/t > 



vi=-=v n =-i [-a + l]i 1/r) 



(37) 



Proof. This follows by making the change of variables U = (1 — Ui)/ui in (35), then writing 

x = (3, y = a. □ 

For A = l k (i.e. 1 repeated k times), k < n, we have that Pi 1 ^^) = efc(z) where efc(z) = 
e k (zi, . . . , z n ) denotes the k-th elementary symmetric function. Noting from (36) that with 

K = l h 

[u]^ = (-l) k (-u;r) k 
we can thus specialize (37) and so reclaim (31). 
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Corollary 4.3 We have 
This is equivalent to (31). 

Proof. It remains to explain the final assertion. This is in fact a consequence of the very 
definition of ifk as given in (12), which gives 

n—k n „ 

<^> # = (-ir fe (rh n (i-^o) 

j=l l=n-k+l 

Substituting (38) in (39) and recalling (10) reclaims (31). □ 

The fact that (31) has been derived independent of (27) means, from the argument of the 
proof of Lemma 4.1, that the difference system (13) can be used to prove (27). 

With (31) substituted in (30), and the entries of A, (14), specialized according to (23), 
we know all terms in (30) except the average — a polynomial in x of degree n/i — which can 
therefore by computed by matrix multiplication. For example, in the case a = (3 = 2, n = 5, 
r = 5, this give for (24) 

23 23x 3197x 2 8993x 3 2117449x 4 793093x 5 



5437500 65250 261000 56550 2035800 203580 
601937x 6 4384a; 7 7457x 8 q in 

+ ^mr-^r + ^>r- f ' x +x (40) 

In general (24) in the case a — (3 must be unchanged (up to a possible sign) by the mapping 
x — 1 — y. One can check that (40) has this invariance. 

We have presented the difference system (12) both for its theoretical interest, and its utility 
in computing the random matrix average (24) in the case /i G Z+. Regarding the latter, the 
case r = 1 of (24), after the change of variables Zi = cos 2 3 -/2, x = cos 2 0/2, and for certain 
values of a, b, corresponds to the /x-th moment of the characteristic polynomial for the classical 
groups Sp(2n), ± (2n) and ± (2n + 1). As such it has appeared in various applied studies in 
random matrix theory [9, 15, 16]. 

It was remarked in the Introduction that the differential-recurrence scheme of [7] can also 
be used to compute the random matrix average (24) in the case \x G Z + . This differential- 
recurrence is equivalent to a (n + 1) x (n + 1) matrix Fuchsian differential equation [11, 20]. The 
fact that there is both a differential and difference system for the Selberg correlation integrals 
is closely related to there being dynamical difference equations associated with solutions of the 
KZ equation [19]; indeed as an example of the latter theory the Selberg integral recurrence 
(10) was reclaimed. In recent years higher order scalar differential equations have been shown 
to characterize certain correlation functions in the two-dimensional Ising model [22, 23, 4] and 
certain generating functions in enumerative combinatorics [13, 14]. Finding characterizations 
of similar problems, outside the class of averages (24), in terms of higher order scalar difference 
equations or matrix difference equations remains an open problem. 
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Another consequence of our results is in providing a fast method of computation of a 
certain class of generalized hypergeometric functions based on Jack polynomials Pi a \z). With 
d a \z) := (a^\K\\/d' K )P!f ! \z) denoting a renormalized Jack polynomial (for the definition of 
the quantity d' K see [8, Eq. (12.37)]), and [u]^ defined by (36), the generalized hypergeometric 
functions of interest are defined by the infinite series 

2 if>K<*2; ^Ei V^fV M. (41) 

In the case n = 1 this reduces to the usual Gauss hypergeometric function. In general, the 
computation of this function from the series is an inherently difficult task due to the need 
to sum over all partitions k [18]. In the special case that a\ is a negative integer the series 
terminates and it is equal to a multivariable polynomial. If furthermore Z\ — • • • — z n — x 
(41) reduces a polynomial of degree n|ai| and it relates to the average (24) according to [8, 
Eq.(13.12)] 



XVx-ZjY) =x N ^ 2 Fi a \-fi,a + (n-l)r;a + /3 + 2(n - l)r; Zl , . . . , z n ) (42) 

, / zi=—=z„=l/x 



Thus the matrix formula (30) can be used to compute this class of 2 F^ using 0(n 3 ) operations. 

In contrast, computation from (41) requires at least 0(e 7r v /2ri / 3 ) operations, due to the sum 
over partitions. 
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A Appendix A — Proof of three-term relations 

Let & n be the symmetric group on {1,2,..., n}, which is generated by the following reflections 
of the coordinates z — (z\, z 2 , ■ ■ ■ , z n 

'■ { z l, ■ ■ ■ z i+l) ■ ■ ■ ) z n) l— ^ { z l, • • • 5 z i+\i z ii ■ ■ ■ 1 z n)i i = 1, 2, . . . , 71 — 1. 

For a function f(z) on C n , we define action of the group & n on f(z) by 

{af){z):=f{a-\z)) for a G & n . 

We say that a function f(z) on C n is symmetric or skew-symmetric if crf(z) = f(z) or crf(z) = 
(sgncr) f(z) for all a G & n , respectively. 

In this section, to specify the number n of variables zi, z 2 , z n , we simply use (p^(z) 
instead of the polynomials (Pij(z) defined in (15). The symbol (n) on the right shoulder of ip it j 
indicates the number of variables of z — (zi, z 2 , z n ) for <fij(z). 
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We define the orbit of the polynomial iplj(z) by 
and define the sum with respect to the orbit by 

which is not always monic. Since, by definition, we have 

(°<Pij) = (<Pij) for any a G 6 

we obtain 



For the point z = (zi, z 2 , . . . , 2 n ) G C n we set 

z% = (zi, . . . , Zi—i, Zi+i, . . . , z n ) G C n 1 i = 1, 2, . . . , n, 

so that, 



¥>S ^(Zn) = (p£j "'(Z!, Z 2 ,..., Zn-i) 



(n-1) 



hJ 



(Zi - Xi) ■ ■ ■ (Zj - Xi) X (x 2 - Zi)--- (x 2 - Z n -i-i) x (x 3 - Z n ^i) ■ ■ ■ (x 3 - Z n _i) 



n— i— 1 



H(z) = H 1 (z) + H 2 (z) 



Lemma A.l Set 

where 

n 

H i( z ) = ~ E ["lO^ _ z k)( x 3 - z k) + <x 2 {x 1 - z k )(x 3 - z k ) + a 3 (xi - z k )(x 2 - z k ) s { ™~ l \z k ), 



k=l 



(43) 



2r 



KKKti 



[xi - 2 fe )(x 2 - 2 fc )(x 3 - z k )s\j 

-(xi - z t )(x 2 - zi){x 3 - 



(44) 



T/ien (H) = 0. 

Remark. Both H\(z) and H 2 (z) are symmetric under the action of (5 n . 

Proof. Put (f)(z) = (f) k (z) := (x\ — z k )(x 2 — z k )(x 3 — z k )s\j 1 \'Zk) in (H) of Lemma 2.1. Since 



(xi - £ fc )(x 3 - z fc ) + (xi - z k )(x 3 - z fc ) + (xi - z k )(x 2 - z k ) sj". 1} (%), 



we have Xlfe=i ^fc 0fc(^) = Hi(z) + H 2 {z) = H(z). On the other hand, from Lemma 2.1, we 
have (ELi V fc = ELi(V fe = 0. Therefore (H) = 0. □ 
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A.l Proof of three-term relation (17) 

In this section we will prove (17) in Lemma 3.3. Throughout this section we assume that 
n < i + j. 

First we will show (Hi) is expanded by (</?l+ij), an d (v^+ij+i)- Since 

ai(x2 - z k )(x 3 - z k ) + a 2 (xi - z k )(x 3 - z k ) + a 3 (xi - z k )(x 2 - z k ) 
= (ai + a 2 )(xi - z k )(x 3 - z k ) + a 3 (xi - z k )(x 2 - z k ) + a x (x 2 - xi)(x 3 - z k ) 

and = (T V 9 !™~ 1 ' > (^fc)! from (43), Hi(z) is written as 

<re6„-i 



= J] J] [("i + a 2)(z k - xi)(x 3 - z k ) + a 3 (z k - xi)(x 2 - z k ) 

k=l <re6„-i 

-ai(x 2 - xi)(x 3 - z k ) a^ j ~ 1 \z k ). (45) 
Under the assumption n < i + j, we have n — i — 1 < j. Then 

(z k - xi)(x 3 - z k )oipf~ X) (%) e 0^l\ d+1 , (z k - Xl )(x 2 - z k )(j^r l \z k ) G 0$ +1 , 

and (x 2 - Xi)(x 3 - 2 fc )o-^ rl) (^fc) e 
Since H\(z) is symmetric, the above implies 

5^ 5^ (^fe - ^i)(^3 - z k )a^r l \z k ) =cs { ^ 



+i,j+i\ z )i 



k=l ae&„-i 



where c is some integer coefficient. In particular, c = 1 is confirmed as follows. The LHS is the 
sum of n x (n — 1)! elements of while - +1 in the RHS is the sum of n\ elements 

of Thus n x (n — 1)! = n! c, i.e. c = 1. In the same way we have 

n 

Yl ( Zk - xi ^ X2 - ^) (7 ^5 _1) ( Sfe ) = s 5+i( z )' 

k=l ci-e6 n _i 

5^ 5^ ( X 2 - ^)(^3 - ^)(7^"" 1) (S fe ) = sj+lj^)- 
fc=l (Tg6„-i 

Therefore, from (45), we expand Hi(z) as 

ifi(z) = (ai + « 2 )s!+i,i+i(^) + "aSiJ+iC^) - "iO^ - (46) 

so that 

(H^ = n\ (ai + afeX^S+ij+i) + a&^SJ+i) - «i(^2 - zi)(</?S+i,j) • (47) 
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Next we will show (H 2 ) is expanded by (</?|+i )t? -) , (fij+i) an d {fl+ij+i)- We first fix k, I 
satisfying 1 < k < I < n. Since, under assumption n < i + j, we have 

PiJ X \ z l7 Z 2 ,..., Z n -i) 

= (zi - Xi) ■ ■ ■ {z n -i-i - Xi) (z n -i - X ± ) ■ ■ ■ (Zj - Xi) 

X (x 2 - Z X ) ■ ■ ■ (x 2 - Zn-i-i) (x 3 - Z n -i) • • • (X 3 - Zj) (x 3 - Z j+1 ) ■ ■ ■ (x 3 - Z n -i), 



n—i—1 



j-(n-i-l) 



n—j—1 



if n < % +j, then for each monomial (Tip\j (%) G O t - ■ one of the following three cases is 
chosen: 

(1) <7</^™ _1 ^(Sfc) has the factor [z\ — X\){x 2 — zi), 

(2) aipfy^^Zk) has the factor {z\ — x 1 )(x 3 — z t ), 

(3) cr^Pi^j l ^{zk) has the factor (x 3 — zi), but does not have the factor [z\ — X\). 

Here we set Si := {a G & n -i ', o satisfies the case (i)} for % = 1,2, 3. Since the numbers of the 
factors (zi-x 1 )(x 2 - zi), (z t - x l )(x 3 - zi) and (x 3 — zi) appearing in aip i ^~ l \z k ) aren-i-1, 
j — (n — i — 1) and n — j — respectively, we have 



#5i = (n 



l)(n-2)!, #5 2 = (i + j-n+l)(n-2)!, #5 3 = (n - j - l)(n - 2)! (48) 



and #Si + + #5*3 = #<5 n _i = (n — 1)!. These numbers are to be used later. 
For each case we will explicitly calculate 

— - — \(x 1 - z k )(x 2 - z k )(x 3 - z k )a(p i ( l ~ 1 \z k ) - (xi - zi)(x 2 - zi)(x 3 - z^aipf 1 '^ (zi) 
z k - zi L 

The point of the calculations below is to use the identities 

1 



z k — Zl 



- z k ) - (a - z,)J = -1 



and 



1 



a - z k )(b - z k ) - (a - z t )(b - zi) = - (b - a) + (a - z k ) + (a- zM. 

Z k — Z[ L J L J 

Case 1 : We assume that cy 9 !™' ^{z k ) is divisible by (x\ — zi)(x 2 — zi). Then we have 

\%) wpt-'Hz,) 



(xi - zi)(x 2 - zi) (xi - z k ){x 2 - z k )' 



Therefore we obtain 
1 



z k - Zi 



(x 1 - z k )(x 2 - z k )(x 3 - z k )o^f A l \z k ) - (x 1 - zi)(x 2 - z t )(x 3 - zi)a^f, 



(xi - z k )(x 2 - z k )cr(Pij 1] {z k ) 



z k — Zi 



(x 3 - z k ) - (x 3 - Zi) 



= (z k - Xi)(x 2 - z k )a(fi\j 1} (z k ) G Og J +1 . 



(49) 



15 



Case 2 : We assume that crip[j ^(z k ) is divisible by (xi — z t )(x 3 — z{). In the same way as Case 
1, we obtain 
1 



Zk ~ Zi 



(x 1 - z k )(x 2 - z k )(x 3 - z k )a^f, l \z k ) - (x 1 - z l )(x 2 - z t )(x 3 - z x )a^f, L \%) 



(n-l), 



h3 



h3 



= (z k - xi)(x 3 - Zk)a(fi\j ^(zk) G OH' 1J+1 . 
Case 3 : We assume that aip^'^^Zk) is divisible only by (x 3 — z t ). Since 



(x 3 - z{) (x 3 - z k ) 



we obtain 
1 



Zk - Zi 



(x 1 - z k ){x 2 - z k )(x 3 - z k )(T(Pij l \z k ) - (xi - zi)(x 2 - zi)(x 3 - zi)aip\™ L> {zi) 



(n-l). 



(x 3 -z k )aip\j l \z k ) 



Zk - Zi 



(xi - z k )(x 2 - z k ) - (xi - zi)(x 2 - Zt) 



-(X 3 - Z k )(T<p!-j l \z k ) \(X 2 - Xi) + (Xi - Z k ) + (x 1 - Zi) 

(z k - x 1 )(x 3 - z k )atp^~ V ' \z k ) + (zi - xi)(x 3 - zi)o^ l ~ x) 



-(x 2 - xi)(x 3 - z k )atp\j ^{zk). 

Note that (x 3 - z^cr^f 1 ^ {%) G and 

fa - z k ){x 3 - Zk^^izk), (Xi - zi){x 3 - zi)a^~ 1] {zi) G 0^} 1J+1 . 
Considering the above cases, since H 2 {z) is symmetric, H 2 {z)/t is expanded as 



(50) 



H 2 {z) 



= c lS \% !j+1 (z) + c 2 s\y +1 (z) - (x 2 - x 1 )c 3 s\l' 1 /z) 



» 



(n) 



(51) 



where ci, c 2 and c 3 are some integer coefficients. On the other hand H 2 (z)/t is also expanded 



as 



H 2 (z) 



l<k<l<n L <tGSi (TGS2 

+ {( Xl ~ Z k)( X 3 - z k)<7<Pi" rl) (%) + { x l ~ z l)( x 3 - 



-(z 2 - a:i)(x 3 - z k )a(p\j l \z k ) 



(52) 



The RHS of (52) includes the sum of 2© (#S 2 +2#S 3 ) elements of 0^\ d+v while s^ 1J+1 (z) in 
the RHS of (51) is the sum of n\ elements of 0^ lj+1 . Thus we have n\ c\ = 2(") (#S 2 + 2#S 3 ). 
In the same way we have n\ c 2 = 2( r J s )#S\ and n\ c 3 = 2(™)#S' 3 . From (48), we obtain 

ci — n — j + i — 1, c 2 = n — i — 1 and c 3 = n — j — 1 . 
Then (51) implies 



(# 2 > = n!r[(n - j + i - 1)(<P™ J+1 ) + {n - i - l)(<p™ +1 ) - (n - j - l)(x 2 - x,)^ 
From (47), (53) and Lemma A.l, i.e., (H 2 ) + (HJ = (H) = 0, we therefore obtain (17). 



(n) 



An) 



.(53) 

□ 
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A. 2 Proof of three-term relation (18) 

In this section we will prove (18) in Lemma 3.3. Throughout this section we assume that 
i + j < n. 

First we will show (Hi) is expanded by (tPij), (^i+ij) anc ^ ( l Pij+i)- Since 

ai(x 2 - z k )(x 3 - z k ) + a 2 (xi - z k )(x 3 - z k ) + a 3 (xi - z k )(x 2 - z k ) 
= -(an + a 2 + a 3 )(z k - xi)(x 2 - z k ) + («i + a 2 )(x 3 - xi)(x 2 - z k ) - a 2 (x 2 - xi)(x 3 - z k ) 

and (43), as we gave the expression (46) of Hi(z) in the previous section, we can expand H ± (z) 
as 

Hi(z) = (ai + a 2 + a 3 )s^ +1 (z) - (oti + a 2 )(x 3 - xi)s^(z) + a 2 (x 2 - Xi)s^ hj (z), 
so that 

(Hi) = n\ (an + a 2 + a 3 )(^!" ) +i) ~ («i + ^)(x 3 - xi)(ip^) + "2(^2 - ziXvS+ij) • (54) 

Next we will show (H 2 ) is expanded by ((fij), (fi+ij) an d (fij+i)- We first fix fc, / satisfying 
1 < k < I < n. Since, under assumption i + j < n, we have 

fij l \ z li z 2i---i Z n -l) 

= (Zx - Xi) ■ ■ ■ (Zj - Xi) 

X (x 2 - Zi) ■ ■ • (X 2 - Zj) (X 2 - Z j+1 ) ■■■(X 2 - Z n _i) (x 3 - Z n ^i) ■ ■ ■ (x 3 - Z n _i), 
N v /s v ,s v ' 

3 (n-j-l)-i i 

if i +j < n, then for each monomial crc/?-" -1 ^^) G 0^~^ one of the following three cases is 
chosen: 

(1) cripf''^ (z k ) has the factor (z x - x 1 )(x 2 - z t ), 



>sr>< 



(2) (T<Pi™j (%) has the factor (x 2 — z{), but does not have the factor (zi — xi), 



(3) c^lj (-Zfc) has the factor (#3 — zj) 

.(*»-!) 

previous section, from (49), we have 
1 



Case 1 : We assume that cr<^"- (2fc) is divisible by (x\ — z\)(x 2 — z{). In the same way as the 



(xi - z k )(x 2 - z k )(x 3 - z k )atp\j l) (z k ) - (xi - zi)(x 2 - z t )(x 3 - zi)atp\j l) (%) 

Z k Z[ l 

= (z k - x 1 )(x 2 - z^a^'V (z k ) e Og } +1 . 



Case 2 : We assume that fy 9 !™' l \z k ) is divisible only by (x 2 — z\). In the same way as the 
previous section, from (50), we have 

— ^ — (xi - z k )(x 2 - z k )(x 3 - z k )a^~ l \z k ) - (xi - zi)(x 2 - z t )(x 3 - z^aip^'^ (%) 

Z k Z\ L J 

= (z k - x l )(x 2 - z k )aip^~ 1] \z k ) + (zi - x l )(x 2 - zi)a^f ~ x) (%) - (x 3 - x l )(x 2 - z k )aip < i n ~ 1) 
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Note that (x 2 — z k )aip\j ^(zk) £ Of* 1 and 

(xi - z k )[x 2 - z k )a^r l \z k ), ( Xl - Zl )( X2 - z^aip^izi) G 0{J +1 . 

Case 3 : We assume that aip^~^ i^z k ) is divisible by (0:3 — zi). Since 

(xi - z k ){x 3 - z k ) = (xi - z k )(x 2 - z k ) - (x 2 - xi)(x 3 - z k ) + (x 3 - x 1 )(x 2 - z k ), 
we have 



(z k - xi)(x 3 - z k )o^f- 1] {z k ) 



h3 



(z k - xi){x 2 - z k )a(p\" ^(zk) + (x 2 - x 1 )(x 3 - z k )<T(Pij L> {z k ) 



(n-1), 



-(x 3 - xi)(x 2 - z k )a^~ 1 \z k ) 

Combining (50) and (55), we obtain 
1 



(55) 



z k — Z\ 



{x x - z k )(x 2 - z k )(x 3 - z k )aip\j l \z k ) - {xx - zi)(x 2 - zi)(x 3 - zi)o^ ^(zj)] 

= (z k - xi)(a;3 - z k )oyf ~ X) (z k ) + (z t - xi)(x 3 - zi)o^ x) {%) 
-(x 2 - Xi)(x 3 - z^aip^f^fa) (from Eq.(50)) 



= (z k - Xi)(x 2 - z k )o^)f- 1] (z k ) - (x 3 - x l )(x 2 - z k )o^f- L) (^ k ) 



(n-1). 



+(zi - xi)(x 2 - zi)a^ 1] (zi) + (x 2 - x l )(x 3 - zi)a(p\j L) {%) 
-(x 3 - x 1 )(x 2 - z^oipf '^Czi) (from Eq.(55)). 



(n-1). 



Here 



(z k - xi)(x 2 - z k )oyf; l \z k ), (z k - xi)(x 2 - zfiatpf) 1] {zi) G {n; 



j'+i' 



(x 2 - ^)cr^ ;j (x 2 - zi)<7(plj \%)eO 

{x 3 -z l )a^- 1 \z l )eOtXr 
Considering the above cases, since the numbers of the factors (zi — Xi)(x 2 — zi), (x 2 — zi) 
and (x 3 — z{) appearing in aip^~ 1 \z k ) are j, (n — j — 1) — i and i, respectively, as we saw in 
the previous section, the function H 2 (z) is expanded in total as 



H 2 {z) = 



2rG)(n-2)!r 



ni 



'jstU z ) + (n-i-j-l) (2s^ +1 (z) - (x 3 - Xl )8^(z)) 
+t(2s^ +l (z) - 2(x 3 - xi)s^{z) + (x 2 - xMlU*)) 



(2n - j - 2)sW +1 (z) - (n + % - j - l)(x 3 - x^s^iz) + i(x 2 - xiK+ij 
so that 

(H 2 ) = n\r[(2n - j - 2)(y^ ) +1 > -(n + i-j- l)(x 3 - Xl )(^) + t(x 2 - n)^- 
From (54), (56) and Lemma A.l, i.e., (H 2 ) + (Hi) = (H) = 0, we therefore obtain (18) 



(56) 
□ 
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B Appendix B — Proof of Corollary 3.4 

We detail the case of Eq. (19) only; the proof of (20) proceeds similarly. 
Induction on j: If j = 1, then Eq. (19) reads 

(an + (k-l)r) (x 2 -x 1 ){(p k>n - k ) = (a 3 + (n- k)r) ((p k -i,n-k+i) + («i + a 2 + (2k - 2)r) ((pk, n -k+i), 

which is confirmed from (17) by setting i = k — 1 and j = n — k. Next we assume the following 
as an inductive hypothesis 

(«i + (k-j + 1)t; r)^^ - ari) J-1 (<?*,,»-*) (57) 

= Hi • )( a 3 + (^-^)^^) i _ i _ 1 (ai + tt2 + (2/c-j>;r).(v?i + fc- i+ i, ri -fe +j -i). 
i=o ^ 2 ' 

On the other hand, from (17) by setting i ^ i + k — j and j — > n — k + j — 1 we have 

(oil + (k- j)r) (x 2 - xi){ip i+ k-j+i,n-k+j-i) (58) 
= (a 3 + (n - % - k + j - l)r) (</? i+fe _ 7> _ fc+J ) + (e*i + a 2 + (i + 2A; - 2j)r) (<p i+fc _ j+1>n _ fc+iJ -) . 

Then the LHS of (19) is written as 

(oil + (k-j)r; t).(x 2 - £i) J (</?fc,n-fc) 
= (ai + (A; - j)r)(x 2 - x x ) 

X H ( • )( a 3 + (^-^)T;r) i _._ 1 (ai + a2 + (2/c-j>;r).(v?i + fc_ i+ i,„-fe +i -i} 

i=0 ^ ' 

(from (57)) 

= H ( • )( a 3 + (n-A;)r;r) i _._ 1 (a 1 + a 2 + (2A;-j>;r). 
i=o ^ 1 ' 

x [(«3 + (n - i - k + j - l)r) (y?; +fe _ iin _ fc+i ) 

+ (ai + a 2 + (i + 2fc - 2j)r) (ip i+k - j+1<n - k+j ) (from (58)) 

= Yl ( ■ ) («3 + (n - A;)r; r) . (ai + a 2 + (2fc - j)r; r) 
i=o V 1 / 

x (a 3 + (n - i - k + j - l)r) ((p i+k _ j:Tl _ k+j ) 

i=i ^ ' 

x (ai + a 2 + (i + 2k - 2j - l)r) (ip i+k _ jin _ k+j ) 

j 

= X («3 + (n - fc)r; r) j _. (a : + a 2 + (2A; - j)r; r) ._ x 
i 

^ J ~ ^ («i + a 2 + (2A; - j + i - l)r) + ^ ~ ^ (ai + a 2 + (i + 2k - 2j - l)r) 

x (<£i+k-j,n-k+j)- (59) 



i=0 
X 
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Here 



J-l 

i 



(ai + a 2 + (2k - j + i - l)r) + _ ^ (an + a 2 + (i + 2k - 2j - l)r) 

J-1 N 



J-l 

i 



+ 



J-l 
i-l 



(cti + a 2 + (2fc - j + i - 1)t) - j i 
y (on +a 2 + (2k-j + i- l)r) - i Q r = Q ( ai + a 2 + (2k -j- l)r) . (60) 
From (59) and (60) we obtain Eq. (19). 
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